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Zitterbewegung is an inherent hallmark of the Dirac equation which arises from the superposition 
between the positive and negative energy part of the spinor states. This exotic quantum phenomenon 
is only observed very recently in traped ions and cold atom systems, but still a challenging task 
for electron systems. Here we present a general theory about the electron Zitterbewegung in a 
topological insulator in the presence of an in-plane electric field. We demonstrate theoretically that 
the topological edge states in topological insulators origin from electron Zitterbewegung and the 
quantum spin Hall effect and quantum anomalous Hall effect can be viewed as a manifestation of 
anomalous electron Zitterbewegung. 

PACS numbers: 03.65. Vf, 73.23.Ad, 72.20.Ee, 72.10.Bg 



The Dirac equation, a combined theory of the quan- 
tum mechanics and the special relativity, predicts the 
existence of the electron spin and the antimatter, which 
were demonstrated experimentally in the last century. 
However, the Dirac equation also predicts some aston- 
ishing effects, for instance, 'Zitterbewegung' and Klein's 
paradox. Zitterbewegung(ZB), an unexpected trembling 
motion of electron, was first predicted by E. Schrodinger 
from the Dirac equation in 1930 The ZB will be com- 
pletely absent if there is no superposition between the 
positive- and negative- energy part of the spinor states. 
This rapid quivering motion around the mean position of 
electron at the speed of light c is proposed to be a pos- 
sible origin of electron spin and magnetic moment 
The ZB is a quantum phenomenon which can even ap- 
pear in non-relativistic cases, such as a particle moving in 
a crystal Q and a Cooper pair in superconductors 
but only obseved experimentally in the trapped ion 
and cold atom 0| systems. However, this prediction is 
still never observed for electrons in free space, because 
this rapid trembling motion displays an extremely high 
frequency loz = 2moc 2 /h ~ 1.5 x 10 21 Hz and a tiny am- 
plitude Xz = h/m$c k, 3.9 x 10~ 4 nm. The detection of 
the oscillation with such high frequency and negligibly 
small amplitude is beyond the reach of the state-of-the- 
art experimental technique. It is natural to imagine that 
the frequency of the ZB should be decreased when the 
energy gap 2moc 2 decreases. 

Semiconductors are good candidates for observing the 
electron ZB phenomenon since its bandgaps range from 
to lOeV, which is much smaller than the energy gap 
(ss lMcV) for electrons in the free space. Very recent 
theoretical works demonstrated that in semiconductor 
quantum structures the electron ZB can be generated 
by the spin-orbit interactions (SOIs) which couples spin- 
up and spin-down subbands and leads to an oscillatory 
motion of the wave packet Q • This ZB behavior is a nat- 
ural consequence of the SOIs which is a relativistic effect 
and generates an entanglement between the spin-up and 
spin-down subbands. The similar jitter behaviors are also 
found in single-wall semiconducting carbon nanotubcs 



and single layer and bilayer graphene [9J . Although there 
are many theoretical works predicting the ZB behavior 
in different semiconductor quantum systems, the exper- 
imental measurement is a very challenging task and ex- 
tremely difficult. This is because that there are many 
electrons (10 10_12 /cm 2 ) in these quantum structure, the 
ZB of each electron shows different phases which smear 
out this oscillation finally. 

The time-reversal invariant topological insulator (TI) 
is a new state of quantum matter possessing an insulat- 
ing bulk and helical metallic edge or surface states, which 
shows a single massless Dirac dispersion at the T point 
of the Brillouin zone. [13, [lit TIs are distinguished from 
a normal band insulator by a nontrivial topological in- 
variant characterizing its band structure. The quantum 
spin Hall effect was proposed independently in graphene 
[l2J and HgTe quantum wells [l3|. The existence of 
edge states and surface states was confirmed by the re- 
cent experiment in HgTe quantum wells 1J| and the 



ang le- resolved photoemission spectroscopy experiments 
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Now we show that the electron quantum dynamics in 
topological insulators gives naturally the origin of the 
edge states based on a general theory of electron ZB. The 
quantum spin Hall effect and quantum anomalous Hall 
effect can be viewed as the manifestation of the electron 
ZB in TIs. 

We consider the single-particle Hamiltonian of electron 
at low-energy regime in the presence of a uniform electric 
field E 



H = e(k)I + d °( k )vi -eE-r, 

i=l 

= e + D° F r, 



(1) 



where <x,-(i = 1,2,3) is the Pauli matrix, e(fc) is the ki- 
netic energy. We assume that the electric field points 
along the x axis E = (E, 0, 0) without loss of any gen- 
erality. The different forms of the off-diagonal term di 
can be used to describe the different important systems: 
1) the two-dimensional electron gas with the Rahsba and 
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Dresselhaus SOIs with d\ — —aky — 0k x , d 2 = ak x +/3k y , 
rf" = 0, where a and (3 describe the strengths of the 
Rahsba and Dresselhaus SOIs, respectively [3[; 2) single 
layer graphene with d\ = vphk x , d 2 = vphk y , d® = 
0; 3) bilayer graphene with d® = -h 2 {k 2 — k 2 )/2m, 
c?2 = —h 2 k x k y /m, d® — [9|; 4) the two-dimensional 
TI HgTc quantum wells with an inverted band structure 
with d\ = Ak x , rig = Ak y , d% = M - Bk 2 (see Figs. 
1(a) and 1(b)) [l3|; 5) three-dimensional TIs, e.g., Bi 2 Sc 3 
and Bi2Te3, whose Hamiltonian would be extended to 
H = Y^,v=a d %( k ) a v® (7 v- eE - r i in wm ch a = J 2X 2, 
and d° m = e(fc), d° 03 = M(k), d° 31 = ^k^Ai = A ^x, 
^12 = Azk y while other d°„s are zero 

mm 

In the following, we show how to derive the electron 
position operator yij(t) evolving with the time t 

y H (t) = e ^ ye -im/H t 

= 2/(0) + (it/h)[H, y] + [H, [H, y}} + • • •, 

= y(0) + (it/h)e« + ^W^T n , (2) 

n—1 

where e* = [e,y] = -id ky e, T x = [D°,y] = 

EiH^d?)^ = = D y , T 2 = [D° 1 Dy} + 

F[D y ,x] = v, ,,/';,/'; Vt,.^: + fd* = /;••••' + i- />■>■. 

• • • , the n-th commutator T n can be obtained from the 
recursive relation, T n — [D°,T n -i] + F[T n _i,x). Gener- 
ally, the analytical expression of the above summation is 
very difficult to obtain because the number of the commu- 
tators increases hierarchically with increasing the order 
n. 

Interestingly, the similarity between the commutation 
relationship [er^er.,-] = 2iEijk&k and the vector cross pro- 
duction (A x B), = EijkAiBjek offers us a new way to 
solve this problem. First, we consider the simple case: 
the electron ZB in the absence of an external electric 
field. In this simple case, the commutation between the 
spin operators is converted into the vector product be- 
tween them. The high-order commutators T n are ac- 
tually the vector product of the low-order T n _i and 
the spinor part Dq of the Hamiltonian H, as shown 
in Fig. 1. We can transfer the algebra summation 
in Eq. ([2]) into the summations of the vector series 
T 2n +i{T 2n )(n = 1, 2, • • • ) utilizing the similarity between 
the commutator [D°, T n ] and the vector product D° x T n , 
T n = [X>o,T„_i] ~ t(D° x T„_i)/2 (see Fig. 1(c)). One 
can see that the vectors T 2n+ i{T 2n ) point along the same 
direction while the lengths of the vectors T n varies with 
the increasing of the order n as T 2n +i — — |-D°/2| 2 T 2 „_i. 
This character allows us to get the analytical expression 
for the electron position operator 

y H (t) = y(0) + j (e y + Ti) + [cos (2\D°\t/h) l] 

+ [sin (2\D°\t/h) - (2\D°\t/h)] , (3) 




FIG. 1: (color online). Schematic of electron Zitterbewegung 
in a HgTe quantum well with an inverted band without and 
with an in-plane electric field (a) and (b). Schematic of the 
recursive relationship between the commutators with different 
orders (c). 



where \D \ = (ElM) 2 ) 1/2 - 

This analytical expression agrees exactly with the pre- 
vious theoretical works for various systems. For a 2DEG 
with the Rashba and Dresselhaus SOIs Q, e y + Ti = 
-ih 2 k x /m - i(aa y - 0a x ), T 2 = {a 2 - f3 2 )k y cr z , T 3 = 
(a 2 - (3 2 )(k x a y - k y a x )T., where £ = a(k x cr x + k y a y ) + 
P{k x a y + k y a x ), \D°\ = (a 2 + (3 2 )k 2 + Aaf3k x k y . For 
a single layer graphene Q, e v + T\ = —ihvpa x , T 2 = 
h 2 VF 2 k y a z , T 3 = h 3 VF 3 k y (k x a y - k y a x ), \D°\ = v F hk. 
The second term e v + T\ represents the classical uni- 
form rectilinear motion, the third and the forth terms 
T 2 and T3 describe the ZB with the oscillating frequency 
2\D°\/h, corresponding to the energy difference between 
the upper and lower energy branches for a given value 
of k. These agreements demonstrate the validity of the 
diagram technique developed by us. 

In the presence of a driving electric field F, the situa- 
tion becomes more complicated. After a lengthy algebra, 
we can also obtain the analytical expression of the ZB 
after more complicated algebra (see the supporting ma- 
terial) 

y H (t) = y(Q) + (it/hy + D^T^V^d, 

n=l 

where the matrix S relates the T n+ \ and T n in the basis 
{D a , D ax P;a / ^C {0,x,y...}}, and Ci is the coef- 
ficient of T\ in this basis. When F =^ 0, S can be di- 
agonalized as S = VEV -1 , and Cq = S _1 Ci, thus the 
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expression of the ZB can be simplified as 

y H (t) = 2/(0) + (it/R)e* - DC 
+ DV exp(iiE/fi)V _1 C . 



400 



(4) 



In general, the theory we developed here can be applied 
to more general cases (see the supporting material) , even 
for the trapped ion and cold atom systems. 

Now we consider electron ZB in two-dimensional topo- 
logical insulators mentioned above including a driving 
electric field. We consider a HgTe quantum well with 
an inverted band structure described by the BHZ model 
13| . The predicted topological edg e states are demon- 
strated experimentally recently [14j . The single-particle 
effective Hamiltonian for the electron is 



H 



D° - Fx, 



(5) 



Bk 2 ) •£, E 



where e = C- Dk 2 , D° = (Ak x ,Ak y ,M 
(a x ,a v ,a z ) T . A, B, C, D, and M are the parameters 
determined by the thickness of the quantum well. In this 
case, D a is expanded in the basis {D°, D x , D y , D x , 
D 0XI , D 0x y, D 0xx \ D xx y, D xxx \ Dy xx2 }, and the S 
matrix in this basis is 
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where Pq.o 


= A 2 k 


2 + 


(M - Bk 2 ) 2 , P , x = - 


2^4 kg 



2iBk x {M-Bk 2 ), P 0ty = -iA 2 k y + 2iBk y (M - Bk 2 ) and 
P _ X 2 = 2B(M - Bk 2 ). 

Next we will show the electron ballistic orbital mo- 
tion for both a topological insulator and a normal band 
insulator, corresponding to the positive and negative pa- 
rameters M. For an electron injected in the Gaussian 
wave packet moving along the x axis, i.e., k x — k$ Xl 
k y = at the beginning t = 0, with spin pointing along 
the z axis, which is perpendicular to the quantum well 



<fi(k x ) 



r i/4 J dk x e 



t) , where T 



V27T 7T 1 

is the width of the wave packet in fc-space and k$ x is 
the initial central wave vector in the x direction for the 
packet. The guiding center of the wave packet (ynit)) 
can be calculated numerically utilizing the explicit ex- 
pression (see Eq. (fj|). The numerical results are plotted 
in Figures 2(a) and 2(b). From these figures, one can 
see clearly the electron orbital motions are very differ- 
ent between the normal band insulator and topological 
insulator phases. For the normal insulator with M > 
(see Fig. 2(b)), the spin-up electron trajectory, i.e., the 
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FIG. 2: (color online). The calculated guiding center (y(t)) 
for a spin-up electron injected in a Gaussian wave packet for 
topological insulator (M < 0) (a) and normal band insulator 
(M > 0) (b). The red and blue curves denote the trajectories 
of spin-up and spin-down electrons, respectively. The relevant 
parameters for the graphene vf and topological insulator A, 
B, C,D, M are taken from Ref. and [20| . An external 
dc electrical field in experiments is applied along the x axis 
E = 1.0 X 10~ 4 mV/nm, i.e., F = 1.0 X 10~ 4 meV/nm, the 
width of the Gaussian wave packet is T 
and the incident momentum is ko x — 0.005nnT 



guiding center (y(t)}, shows an oscillation against the 
propagating time as in a conventional 2DEG with SOIs; 
while for the topological insulator M < (see Fig. 2(a)), 
a spin-up electron shows a surprising behavior trajectory 
which bends opposite to the y axis monotonously. This 
behavior make it possible to observe the Zitterbewegung 
first time in a solid. As we mentioned before, electrons in 
the normal band insulator oscillate with different phases, 
leading to smearing out of the ZB. 

' Now we turn to demonstrate theoretically the edge 
states in a spin Hall bar with a finite width. For 
an realistic case, i.e., the HgTe quantum well sample 
with finite size, one can image that the electron will be 
bounced back when it reaches the boundary of the spin 
Hall bar. From the expression of the spin-orbit force 
f v ~[H,[H,y]]~A[M + B{k 2 y - k 2 )}a x - 2ABk x k v a v - 
A 2 k x a z , one can see that (the last term — A 2 k x a z ) the 
spin-orbit force always bends the spin-up (down) electron 
downward (upward). Therefore the electron bounced 
back from the boundary will bend to the boundary again, 
and forms a edge channel in the vicinity of the boundary 
finally. This spin-orbit force could make the edge skip 
state robust against the non-magnetic impurity scatter- 
ing. The spin-orbit force will lead to the spin Hall effect 
in a 2DEG in the presence of SOIs, Notice that the en- 
ergy of incident electron here is (E) = C + M — (B + 
D)(T 2 + k^ x )/2, which locates inside the bulk gap. 

To explicitly get the edge states in quasi-one- 
dimensional systems, we can map the infinite 2D system 
(x, Y G (— oo, +oo)) into a quasi-one-dimensional system 
with a finite width (x,y G (— 1,+1)) using the trans- 
formation ±l/(l=py) = Y ± 1, where y and Y are the 
dimensionless coordinates for the infinite and finite sys- 
tems, respectively. For Y > (for Y < case, we can 
adopt the same transformation — 1/(1 + y) = Y — 1 ) 
ky = —idy = (y — l) 2 k y , then the original 2D Hamilto- 
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FIG. 3: (color online). The trajectories of the electron Gaus- 
sian wave packets in a spin Hall bar mapped from 2D infinite 
system. The inset shows the transverse confining potential. 
For the topological insulator (M < 0), the edge state appears 
in the vicinity of the QSH bar; while for the normal band 
insulator (M > 0), the trajectory shows a quivering behavior. 



nian can be written as 



H% = Ak x a x + Ak y a y + (M - Bk 2 x - Bk\)o z , 
= Ak x a x + A'k y a y + (M - Bk 2 x - Bk 2 )a z 



-B{{y-l) A + 2{y-lf-l}ky 

V(y), 



H'y 



(6) 



where A' — A(y — l) 2 and the last term of Eq. [BJ 
V(y) = ~B{{y - l) 4 + 2(y - l) 3 - l}k 2 y a z = -Blty*v{v) 
can be viewed as the confining potential without any 
loss of the generality. Thus a spin-up electron will feel 
a confining potential V(y), and the ZB solution can be 
obtained by converting the 2D solution into this con- 
fined system utilizing the above coordinate transforma- 
tion. The spin-up and -down edge states appear at the 
vicinity of the opposite boundaries of the HgTe spin Hall 
bar in the topological phase (M < 0) when y — >■ ±1. For 
the normal band insulator phase (M > 0), one can see 
that the electron exhibits an normal oscillation in the bar 
(see Fig. 3). 

Finally, we would like to point out that the quantum 
anomalous spin Hall effect can also be understood from 
the electron ZB in TIs. Considering a ferromagnetic TI, a 
giant Zecman splitting results in a spin-up inverted and a 
spin-down normal band structure, we assume these mag- 
netic ions polarized along z axis. This situation corre- 
sponds to the BHZ Hamiltonian with the positive and 
negative M for spin-up and spin-down electrons, respec- 
tively. From the above discussions, the spin-up (down) 
electrons show normal oscillating (bending) trajectory 
(see Fig. 2). Therefore, the edge states only appear for 



spin-down electrons, i.e., the quantum anomalous spin 
Hall effect. Similar behavior can also be found in three- 
dimensional topological insulators, because the form of 
the Hamiltonian of 3D TI is almost the same as that of 
2D TI. Therefore one can expect that the bending tra- 
jectory of electrons in 3D TIs will lead to the topological 
surface states. 
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Calculation for yn{t) in spin-orbit coupled systems 

We consider the single-particle Hamiltonian of electron in the presence of a uniform electric field E 

H = e(k)I + d°{k)a l -eE-r 

= e + D°~F-r, (SI) 

where <jj(i = 1, 2, 3) is the Pauli matrix, e(k) is the kinetic energy. We assume that the electric field points along the 
x axis E — (E, 0, 0) without loss of any generality. 

In the following, we will derive the electron position operator yii{t) evolving with time in the Heisenberg picture. 

yH (t) = e ^t/H ye -im/n 

= 2,(0) + iit/h) [H, y] + [H, [H, y]] + ••• 

= ymm y + ^^ + ^^T n , (S2) 



n=l 



here we show the notions explicitly: 

= [eV,x]=d Kky e, (S3) 

T 1 = [D\y] = {-id ky dP l )^ 

= d v i a i =D y , 
T 2 = [D ,Dy}+F[Dy,x} 

= d°d y {v\<T J }+F(-id k Ji)a l 

= 2/ , (F'lf'rr'' + Fd^G 1 

= D 0xy + FD yx , ... (S4) 
The n-th commutator T n can be obtained from the T" n _i, 

T n = [D°,T n _ 1 ] + F[T n _ 1 ,x], (S5) 

and: 

[D°,D a ] = D 0xa , 

[D axb ,x] = D axxb + D axbx , (S6) 
£,ax(f>xc) = (2i) 2 (P a , c D b - P a , b D c ) , 

where a, b, c C {0, x, y, x 2 , xy, • • • }, which can be any combination of x, y; P a .b = D a ■ D b = dfd b is the "dot" product 
of D a and D b . 

The above relations indicate that all T n are linear combination of {D a , D axb ; a ^ b C {0, x, y, x 2 , xy, ■ ■ ■ }}, the 
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recursion relation between T n and T„_i can be represented by the matrix S in the vector space spanned by {D a , D bxc }: 

T n = C n ^ a [D°,D a ]+FC n ^. a [D a ,x} 

a—a,bx c 

= C n -i ja D ()xa + FC n -i ja D a ' x 



a—a,bx c 



or in the compact matrix form: 



= E E S Pa C n - lia ) DP, (S7) 

0—a,bxc \a—a,bxc 



T n — DC n — D-SC„_i 

= ••• = D-S"- 1 Ci, (S8) 

r> = (D°,D*,Dy,.-.), 

C n — (Cn,0? Cn.xi Cn,y? * ) • 



thus 



y H (t) = y(0) + (it/ny + ^4 



2 

+ T) y W^T s n-l Ci ^ (sg) 

When F = 0, i.e., in the absence of an electric field, can be calculated without diagonalize S because 

S 2fe + 2 = P ,oS 2fc = P fe S 2 , 

so that 

(itlt)\ 2 f x y 

y H (t)=y(0) + ( l t/h)ey+ ( -^ 

sm(itPn(?/h) „ 
+ °- 0/ y D-S 2 C 1 

M),0 

+ ( p 3/^ - -ijD-SCx, (S10) 

this result is consistent with previous results with zero field. 
When F ^ 0, S can be diagonalize as 

S = VEV"\ 

and 

y H (t) = y{0) + (it/hW 

-D-Co + D- Vcxp(iiE/ft)V _1 Co, (Sll) 

where C = S^Ci. 

Generalization of the calculation 

Consider the generalized Hamiltonian under an external force field T r 

H = e + V°- J" r , (S12) 
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where 2?° is the generalized spin-dependent operator in k space, while belongs to the closed algebra A satisfying: 
for any V a ,V b 6 A 

[V a ,V h ]=C abc V c , (S13(a)) 
[V a ,T r ]=B ab V\ (S13(b)) 

where C a bc is the structure constant (/c-dependent) of the closed algebra A; B ab is the structure constant of the 
extended closed algebra include both A and T T . 
In the same way, 

T n = D-C n = S3-SC„_i 

= ••• = S>-S n_1 Ci ) (S14) 
S = (Z>°,2? X ,X>V--), 

(Cn.O: ^n,x: Cn.yi ' ) ■ 

and 

y H (t) = 2/(0) + (it/fi)e» + ^ 

+ M^S^d. (S15) 

^-^ TV. 

n=l 

In spite of the position operator, other operators like velocity and spin can be calculated in the similar way. 



